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C/3 ■ ABSTRACT 

(-h ' Wc present Li and Fc abundances for 87 stars in the globular cluster M4. obtained by 

O I using high-resolution spectra collected with GIRAFFE@VLT. The targets range from 

1 ■ the Turn- Off up to the Red Giant Branch Bump. The Li abundance in the Turn-Off 

stars is uniform, with an average value equal to A(Li)= 2.30±0.02 dex (er= 0.10 dex), 
consistent with the upper envelope of Li content measured in other globular clusters 
and in the Halo field stars, confirming also for M4 the discrepancy with the primor- 
dial Li abundance predicted by WMAP+BBNS. The global behaviour of A(Li) as a 
function of the effective temperature allows to identify the 2 main drops in the Li evo- 
' lution due to the First Dredge-Up and to the extra- mixing episode after the Red Giant 

Branch Bump. The measured iron content of M4 results to be [Fe/H]= -1.10+0.01 
dex (<t= 0.07 dex), with no systematic offsets between dwarf and giant stars. 
The behaviour of the Li and Fe abundance along the entire evolutionary path is in- 
compatible with theoretical models including pure atomic diffusion, pointing out that 
an additional turbulent mixing below the convective region needs to be taken into 
account, able to inhibit the atomic diffusion. The measured value of A(Li) and its ho- 
mogeneity in the Turn-Off stars allow to put strong constraints on the shape of the Li 
profile inside the M4 Turn-Off stars. The global behaviour of A(Li) with the effective 
temperature can be reproduced with different pristine Li abundances, depending on 
the kind of adopted turbulent mixing. One cannot reproduce the global trend start- 
ing from the WMAP +BBNS A(Li) and adopting the turbulent m i xing d escribed by 
I^S , iRichard et alj (|2005| ) with the same efficiency used bv lKorn et al.l (|2006af ) to explain 

■ the Li content in NGC 6397. In fact such a solution is not able to well reproduce 

simultaneously the Li abundance observed in Turn-Off and Red Giant Branch stars. 
Otherwise, the WMAP+BBNS A(Li) can be reproduced assuming a more efficient 
turbulent mixing able to reach deeper stellar regions where the Li is burned. 
We conclude that the cosmological Li discrepancy cannot be easily solved with the 
present, poor understanding of the turbulence in the stellar interiors and a future 
effort to well understand the true nature of this non-canonical process is needed. 

Key words: stars: abundances - Stars: atmospheres - Stars: Population II - 
(Galaxy:) globular clusters: individual (M4) 



* Based on observations collected at the ESO-VLT under pro- 
gram 081.D-0356. Also based on observations collected at the 
ESO-MPI Telescope under program 69.D-0582, and with the 
NASA/ESA HST, obtained at the Space Telescope Science In- 
stitute, which is operated by AURA, Inc., under NASA contract 
NAS5-26555. 



1 INTRODUCTION 

The determination of the initial Li abundance in low-mass, 
metal-poor, Population II stars is one of the most debated 
astrophysical topics. Li is one of the few elements syn- 
thesized during the primordial Big Bang nucleosynthesis 
(BBNS), and its initial abundance in metal poor stars sets 
constraints to the BBNS. 
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Li is a very fragile element, that is destroyed when 
the temperature is larger than ~2.5 10 6 K through the 
7 Li(p,a) 4 He reaction. Hence, whenever Li is produced 
in the stellar interiors during hydrostatic burnings, it is 
immediately disintegrated. The discovery of a constant 
Li abundance in une volved, Population II stars, the so- 
called Spite Plateau (|Spite fc Spitel 1 1982T ) - has been in- 
terpreted as the signature of the primordial abundance 
of Li, produced during the BBNS. When the lithium 
abundance is expressed as A(Li)=log(n(Li)/n(H)) + 12, the 
Spite Plateau turns out to be betwe en 2.1 and 2.4, de- 
pend i ng on the adopted T e t f scal e dBonifacio fc Molard 
19971; ICharbonnel fc Primas i [20051; lAsplund et all 1200^ : 



dwarf stars have been until now limited to t he few closest 
objec ts, namely 47 Tuc dBonifacio et al. 20071; D'Orazi et al.l 
2010T). NGC 6397 dBonifacio et all 12003; iKom et alj|2006al ; 
Gonzalez Hernandez et alj 120091 ; iLind et alj l2009bh . NGC 



Bonifacio et al.ll2007l ; lAoki et alJl200gf T 

The recent estimate of the cosmological baryon density 
obtained with the Wilkinson Microwa ve Anisotropy Probe 
(WMAP) satellite dSpergel et alj2007h from the power spec- 
trum of the cosmic microwave background fluctuations, 
throws into a crisis the classical interpretation of the Spite 
Plateau. In fact, by combining the WMAP results with stan- 
dard BBNS calculations, the most recent esti mate of Li is a 
prim ordial abundance A(Li)= 2.72±0.06 dex (|Cvburt et al.l 
[2ooi)FI. This value is significantly higher, by at least a factor 
of 3, than the Li abundance derived from metal-poor dwarfs. 

At present, the WMAP+BBNS / Spite Plateau discrep- 
ancy is still unsolved and different possible solutions have 
been advanced and investigated: (1) an inadequacy of the 
standard BBNS treatment; (2) Population III stars could 
have been capable to destroy some of the pristine Li 
and Population I I stars were born from a Li-depleted gas 
dPiau et al.|[200f3 ) ; (3) Li may be depleted in the photo- 
spheres of Population II stars (born with WMAP+BBNS 
Li abundance) by atomic diffusion: this process alters the 
surface chemical composition of the star, hiding a fraction 
of Li (but also He, Fe and other heavy elements) below the 
observable photospheric layers. 

Galactic globular clusters (GCs) provide an alternative 
way to investigate the primordial Li abundance, for they 
host samples of stars approximately coeval, with uniform ini- 
tial abundances of several elements - Fe among them - and 
whose evolutionary status is easily determined from pho- 
tometry. Also, GCs offer the possibility to trace the diffu- 
sion effects along the entire evolutionary path of Population 
II stars. In fact, the amount of diffusion can be easily mea- 
sured from the offset of the surface abundance of elements 
like Fe - that are not involved in nuclear burnings - between 
the turn off (TO) stars (where the diffusion affects apprecia- 
bly the photospheric abundances) and the red giants (where 
convection restores the original photospheric abundances). 
Systematic differences between the Fe content of Sub Gi- 
ant Branch (SGB) and red giant branch (RGB) stars have 
been detected in M 92 dKing et all 1 19981 ) and in NGC 6397 
dLind et alj|2009bl ). On the other hand independent studies 
have found a high level of homogeneity in Fe an d other ele- 
ments in th e latter cluster dGratton et al. 200 and in oth - 
ers (see also lCastilho et al.ll20obl ; ICohen fc Melendez f2005l ). 

Extensive chemical screenings of Li content in GCs 



1 Ado pting different rates for the 3 He(o,7) 7 Li reaction. ISteigmanl 
ll2007f ) derived a slightly lower value, A(Li)=2.65, still in disagree- 
ment with the Spite Plateau. 



6752 dPasauini et alj|2005l ) and M92 dBoesgaard et al.|[l99gl ; 
lBonifacidl2002r i. The metal-poor GC NGC 6397, in partic- 
ular, has been recently the object of an intense debate re- 
garding the discrepancy betwee n co smological and stell ar Li 
abundance. iKorn et alj (|2006al ) and lLind et~aH d2009ri ) find 
a trend between A(Li) and T e / / along t he SGB and the early 
RGB, roughly consistent with models (|Richard et al.l [20051 . 
hereafter R05) that start with the BBNS A(Li) and include 
atomic diffusion moderated by some ad-hoc turbulent mix- 
ing that brings back to the surface some of the Li (and other 
elements) diffused below the c onvective envelope during 
the m ain sequence (MS) phase. iGonzalez Hernandez et ahl 
(|2009f ) have studied a sample of MS and SGB stars in the 
same cluster, and found A(Li) decreasing with decreasing 
T e ff in both groups of stars, with the SGB stars showing 
systematically higher values of A(Li) at a given T e //. This 
pattern cannot be reproduced satisfactorily by any existing 
theoretical model starting with the BBNS value of A(Li). 

It is relevant to highlight that the investigation of A(Li) 
in GCs needs some special care, because of the occurrence 
of chemical anomalies in the light elements. In fact, all the 
GCs studied so far exhibit large star-to-star variations of 
C, N, O, Na, Mg and Al and in particular a nt i-correlations 
betwe en O and Na, and between Mg and Al dCarretta et al.l 
l2009al ibh. The widely accepted scenario to explain these pat- 
terns is that a second stellar generation was born in the 
cluster from the hydrogen-processed material ejected by the 
first population jGratton et al.ll200ll ). In this framework, Li 
could have a relevant role. A crucial point to recall is that 
the NeNa cycle, responsible for the pollution of the cluster 
gas, occurs at temperatures of 4-5T0 7 K, 20-30 times higher 
than the Li-burning temperature. Thus, unless the pollut- 
ing stars are also able to produce Li, the second generation 
stars should be Li-poor with respect to the first one and 
Li-O correlations and Li-Na ant icorrelat ions are expected. 
Previous in vestigations about th is issue lead to conflicting 
conclusions: IPasauini et all (120051 ) observed for the first time 
the Li-Na anticorrelation (and also the Li-O correlation) in 
9 TO stars o f NGC 6752 and a sim ilar behaviour has been 
found also bv lBonifacio et al.l 120071 ) in 4 TO stars of 47 Tuc. 
In their sample of NGC 6397 TO stars, ILind et alj (|2009bl ) 
detected 3 stars out of 100 with low Li and high Na abun- 
dances, suggesting a possible Li-Na anticorrelation, while 
all the other TO stars share the same Li abundance over a 
ran ge of ~0.8 dex in A ( Na). R ecently, iD'Orazi et all d,2010l ) 
and lD'Orazi fc Marind (|2010l . hereafter DM10) presented Li 
abundances in TO stars of 47 Tuc and RGB stars of M4, re- 
spectively. In both cases, first and second generation stars 
share the same Li abundance, regardless of the Na and O 
content, but with different star-to-star scatter between the 
two stellar groups. 

In this paper, we present the Li abundance measured 
from high-resolution spectra of a sample of stars ranging 
from the TO to the RGB Bump in the GC M4. This is 
the first determination of Li content in the dwarf stars of 
this cluster. The Li content in the TO stars of M4 is then 
compared with the Li abundances measured in other GCs, 
in the field stars, and with the WMAP+BBNS predictions. 
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The behaviour of A(Li) as a function of [Fe/H] is discussed 
in the light of stellar evolution models with and without 
atomic diffusion. Also, the Fe content along the entire evo- 
lutionary path is discussed and it is used to quantify the 
efficiency of atomic diffusion. Finally, we discuss the possi- 
ble role of turbulent mixing in the theoretical models and the 
open questions related to the determination of the pristine 
Li abundance of M4. 



2 OBSERVATIONAL DATA 

This work is based on high-resolution spectroscopic ob- 
servations of 87 stars in M4 performed with the GI- 
RAFFE@FLAMES spectrograph mounted at VLTQESO, 
within a project aimed to sea rch for anomalous a bundances 
in the BSS atmospheres (see lFerraro et all [20061. All stars 
have been observed with the HR15N, HR18 and HR22 grat- 
ings that include the Li doublet at 6707.8 A, the H a Balmer 
line and several Fe I lines. The spectral sample includes 
5f TO/SGB stars and 36 giants up to the RGB Bump 
magnitude level. Target selection, reduction of the spectra 
and measur e ment of the radial velocities are discussed in 
ILovisi et al.l (|2010T ). The position of the targets in the (V, 
V-I) colour-Magnitude Diagram (CMD) is shown in Fig. [1] 
(left panel). The Signal-to-Noise Ratio (SNR) around the 
Li line is ~50-80 for dwarf stars and rises up to about 300 
for the brightest RGB stars. The membership of all the tar- 
gets has been confirmed by the combined information of 
proper motions, radial velo cities and metallicities, as dis- 
cussed in ILovisi et all l|201Ch . Only one target in our sample 
ha s been found to b e in co mmon with the list of binaries 
by ISommariva et all (|2009l ) and it was excluded from the 
analysis. 



3 ANALYSIS 

3.1 Atmospherical parameters for TO/SGB stars 

Atmospherical par ameters for SGB stars have been al- 
ready presented in ILovisi et al.l (|2010h . Here we briefly dis- 
cuss in more details only the derivation of T e //- Tem- 
peratures have been computed from a \ 2 minimization 
between the observed wings of the H a Balmer line and 
a grid of synthetic spectra. Synthetic spectra have been 
computed by means of the SYNTHE code from a grid 
of ID LTE plane-parallel model atmos pheres, obtain ed 
with ATLAS9 code in its Linux version (|Kurucj[l993al lb1; 
ISbordone et al.l I2004F ). Adopted model atmospheres are 
based on the NEW Opacity Dis tribution Functions de- 
scribed bv lCastelli fc Kurucj (|2003r i without inclusion of the 
approximate overshooting and with a-enhancement abun- 
dance patterns of [a/Fe] = +0.4 dex. For the modeling of 
the Balm er line, we adopted the se lf-resonance broadening 
theory bv lAli fc Grieml (|l965l , I1966T ). 

The dominant source of error for T e f / derived from the 
Balmer line is due to the residuals in flat-fielding and blaze 
function correction, that can introduce residual and spu- 
rious curvatures in the spectrum, afflicting the continuum 
placement. Such an effect is not easy to model, but in the 
fiber-fed GIRAFFE spectra it is less critical than in the slit 




V-I A (Li) NLTE 

Figure 1. Left panel: colour-magnitude-diagram of M4 (gray 
dots) with spectroscopic targets marked as black filled circles. 
The solid line is the best-fit isochrone (see text for details) . Right 
panel: behaviour of A (Li) as a function of luminosity. The two 
dashed horizontal lines denote the position of the two detected 
drops of A(Li) (see text for details). 



spectra, since the flat-field is acquired through the same op- 
tical path used for the stellar spectra. In order to minimize 
this uncertainty, observed and synthetic spectra have been 
normalized by using the same continuum regions, typically 
±f00 A from the line centre. 

In order to assess the internal error in the T e //, we re- 
sorted to Monte Car lo simulations, similar t o the procedures 
alread y adopted by iBonifacio et al.l (|2007T ) and iLind et all 
(|200ST) . The effect of the photon noise is quite small: we per- 
formed Monte Carlo simulations of 1000 synthetic spectra 
around the H a line injecting Poissonian noise in order to 
reproduce the typical SNR of our spectra and we repeated 
the procedure to derive T e //. We find that at SNR= 50 the 
error is of about 20 K, due to the large number of pixels used 
in the x 2 minimization. A similar check has been performed 
to assess the impact of the adopted fitting windows. Finally, 
we attribute to our T e // random errors of the order of 30-40 
K. Other sources of uncertainties can affect the Balmer T e // 
and considered as systematic errors, for instance the gravity 
dependence of H a , the residual curvature in the flat-fielding, 
the adopted bro adening theory and f inite pixel-step in our 
spectra (see e.g. IBonifacio et al. I l2007l . for a detailed discus- 
sionof thes e source of error s ). In p articular, IBonifacio et al.l 
|2007h and lSbordone et al.l (|201dl ) discussed the sensitivity 
of the H a line to the gravity and found that it is important in 
the low metallicity regime ([Fe/H]<-2.0 dex), significantly 
lower than that of M4. For the TO stars we estimate a vari- 
ation of only ~10-20 K in T e // for a variation of 0.5 in log 
g, while for the coldest SGB stars such a sensitivity rises to 
~40 K. Bearing in mind these sources of uncertainty and we 
conservatively attribute a typical, absolute error of 100 K to 
TO/SGB stars and of 150 K to the coldest SGB stars. 
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3.2 Atmospherical parameters for RGB stars 

For RGB stars T e ff cannot be inferred with the proce- 
dure described above, since, the H a wings are not sensitive 
enough to T" e // variations in cold stars. Some authors use 
the wings of H a line to derive T e ff also for RGB stars, but 
it is worth noticing that this method, when applied to giant 
stars, is subject to very high uncertainties (of 300 K at least). 
For stars with V<15.3, T e // have been derived by project- 
ing their position in the CMD on the best-fit isochrone. (left 
panel of Fig. [T] where the CMD is corrected for differential 
reddenin g). On the basis of the standard procedure (see for 
example Piotto et al.lll999l ). we estimate that the differen- 
tial reddening is of the order of <SE(B-V)= 0.25 in this clus- 
ter. The isochrone with appropriate metallicity (Z= 0.004, 
corresponding to [Fe/H] = — 1.01 for [a/Fe]=+0.4) and age 
(12 Gyr), has been sel ected from the BaSTI database □ 
(|Pietrinferni et al]|2006h , and transformed into the observed 
plane by adopting E(B-V) = 0.32 and (m-M) g= 11.30, con- 
sistent with the results by iBedin et al.l (|2009l ) . It is worth 
noticing that because of the high value of the extinction, 
we applied a reddening correction that takes into account 
the dependence of A\/Av on the star surface gravity, T e // 
and [Fe/H T. We did this by fol lowing the same procedure de- 
sc ribed in IBedin et al.l (l2009l l. adopting the extinction law 
of lCardelli et al. I (| 19891 ). re-scaled to R y= 3.8, which is the 
value measured in the direction of M4 l|Hansen et al.l [20041 ; 
IBedin et al]|2009h . 

The colour- T e / / transformations of the BaSTI database 
have been computed from the same model atmospheres 
adopted in our spectroscopic analysis. Hence such a T e ff 
scale is formally homogeneous with that adopted for the 
dwarf stars. Indeed the BaSTI isochrone T e ff scale well 
agrees with that derived from the H a line for the TO/SGB 



stars: we find an average difference T e ^ ST/ - 



T 



eff 



-28 K 



(o"=126 K) for stars with V>15.3. Thus, no relevant system- 
atic offset is introduced between the abundances computed 
for dwarf and giant stars (the effect of different T e ff scales 
on the derived abundances is discussed in Sect. I6.2|l . Tak- 
ing into account the photometric error, the uncertainty in 
E(B-V) and that arising from the differential reddening cor- 
rection, we attribute a typical error of 100 K to T e / / of RGB 

stars. 

Similarly to the dwarf stars (see lLovisi et al.ll2010t ). the 
gravities have been obtained from the position of the stars in 
the CMD, and the microturbulent velocities have been de- 
rived spectroscopically, by erasing any trend between A(Fe) 
and the equivalent width (EW). Main parameters, Li and 
Fe abundances for all the target are listed in Table 1, the 
complete version being available in electronic form. 

3.3 Lithium 

Li abundances have been measured by comparing the ob- 
served EWs of the Li I resonance doublet at 6707.8 A with 
the EW computed from a grid of synthetic spectra. A curve 



of growth has been built for each star and we derived A(Li) 
interpolating each curve at the observed EW value. Due to 
the small wavelength separation of the Li components, this 
line can be reasonably approximated with a single gaussian 
profile at the GIRAFFE resolution. EWs have been mea- 
sured with our own FORTRAN procedure, that performs 
a gaussian fit on a spectral window selected interactively. 
The continuum level is traced locally by considering a re- 
gion of ~20 A around the Li line, rejecting spikes, cosmic 
rays and spectral features with a a-clipping algorithm and 
adopting as continuum level the peak of the flux distribu- 
tion computed from the surviving points. As sanity check we 
performed an independent EW measurement by using the 
IRAF task splot, finding an average difference EWiraf- 
EWfu— -0.09 mA with a dispersion a— 2.35 mi. Such a 
small difference is irrelevant to the purposes of this work. 

Atomic data for the Li doublet, including its hyperfine 
structur e and isotopic splittin g (with N( 6 Li)/N( 7 Li)= 8%), 
are from lYan fe Drake! (|l995h . Corrections fo r NLTE effects 
have b een applied by interpolating the grid of lCarlsson et all 
(|l994h : derived NLTE corrections are quite small, of the or- 
der of 0.01 dexQ. 

The error in A(Li) is computed by adding in quadra- 
ture the two main sources of uncertainty, namely the errors 
arising from T e // an d from the EW measurement. A (Li) 
turns out to be very sensitive to the used T e //, because Li 
is almost totally ionized. The typical error in T e // (±100 K) 
translates in an error in A(Li) of ±0.07 dex for TO/SGB 
stars and slightly higher (±0.09 dex) for RGB stars. It is 
customary to estimate th e error in EW measurement by 
using the Cayrel formula (|Cavrellll988h . This gives ±4-4.5 
mA for TO/SGB stars and ±1-2 mi for RGB stars, trans- 
lating in A(Li) errors of the order of ±0.05-0.07 dex and 
of ±0.01-0.02 dex, respectively. The Cayrel formula neglects 
the uncertainty due to the continuum location, that can have 
a relevant impact in spectra with low SNR. We estimated 
the continuum location uncertainty from the average spread 
of the residual spectrum around the Li line and performing 
several tests with different assumptions of the continuum 
level. Finally we derive an uncertainty of ±0.10 dex for the 
TO/SGB stars with lowest SNR and a few hundredths of dex 
for the RGB stars. Other sources of errors (linked to gravity 
and microturbulent velocity) can be assumed to be negligi- 
ble since they correspond to an error in A(Li) of ±0.01 dex 
each. Final absolute errors in A(Li) range between ±0.08 
and ±0.15 dex, according to the stellar evolutionary stage. 



3.4 Iron 

Our spectral dataset allows to measure the Fe abundance in 
both TO/SGB and RGB stars. In turn these measures can 
be used as tracers of the occurrence of atomic diffusion pro- 
cesses along the evolutionary path. Fe abundances for all the 
target stars have been derived from the EW measurement by 
using WIDTH code in its Linux version (ISbordone et al.l 



2 http:/ /albione. oa-teramo.inaf.it/ 

3 The adopted isochrone is calculated without atomic diffusion. 
Including atomic diffusion would change the age of the isochrone 
of about 1 Gyr but the T e y y-colour relation would be essentially 
unchanged. 



4 Note that the adoption of the NLTE corrections of lLind et al.l 
ll2009ah would have a negligible (of about a few hundredths) im- 
pact on the Li abundances. 

5 We have employed a version of the original code by 
R. L. Kurucz modified by F. Castelli, in order to use 
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Table 1. Identification numbers, coordinates, temperatures, gravities, EW of the Li line, 
A(Li) and [Fe/H] abundances, SNR around the Li line for the observed stars. A complete 
version of the table is available in electronic form. 



ID 




RA 




Dec 


T eff 


log g 


EW(Li) 


A(Li) 


[Fe/H] 


SNR 




(J2000) 


(J2000) 


1 \,"\ 

(K) 


(dex) 


(m.A) 


(dex) 


(dex) 




8460 


245 


9471429 


-26 


3749119 


5000 


2.8 


22.7 





96 


-1.14 


104 


8777 


245 


9373405 


-26 


3613910 


5150 


3.2 


11.8 





79 


-1.13 


85 


9156 


245 


9069486 


-26 


3438488 


4950 


2.7 


20.9 





85 


-1.06 


217 


13282 


245 


7901731 


-26 


3909501 


5850 


3.9 


39.5 


2 


11 


-1.00 


53 


28007 


245 


8049983 


-26 


6687472 


5100 


3.0 


13.6 





81 


-1.16 


116 


28890 


245 


8083640 


-26 


6339178 


5850 


3.8 


30.5 


1 


98 


-1.06 


62 


28999 


245 


7440673 


-26 


6298102 


6000 


1.0 


34.4 


2 


16 


-1.20 


68 


29074 


245 


8010945 


-26 


6273020 


6150 


3.9 


34.8 


2 


28 


-1.00 


18 


29411 


245 


7990049 


-26 


6149585 


5500 


3.7 


15.9 


1 


34 


-1.05 


60 


29643 


245 


7572145 


-26 


6058864 


5800 


3.8 


40.3 


2 


08 


-1.19 


68 


29729 


245 


7547183 


-26 


6022333 


5950 


3.9 


46.8 


2 


28 


-1.18 


58 


30167 


245 


7958060 


-26 


5866438 


5950 


3.8 


35.6 


2 


14 


-1.04 


60 


30475 


245 


7195282 


-26 


5759641 


5100 


3.0 


12.0 





71 


-1.05 


158 


30922 


245 


7432036 


-26 


5602858 


6050 


1.0 


33.1 


2 


18 


-1.16 


48 



|2004 ). We employed two different linelists for TO/SGB and 
RGB stars, in order to include only unblended transitions, 
by taking into account the different degree of blending and 
line strength between dwarf and giant stars. The number of 
Fe I lines is limited by the relatively small wavelength cover- 
age of our spectra (we remark that the spectral regions sam- 
pled by the gratings HR18 and HR22 are strongly affected 
by deep telluric absorptions). Finally, we derived the abun- 
dance from ~6-14 and -18-22 Fe I lines, for TO/SGB and 
RGB stars respecti vely. Atomic data are from th e classical 
line compilation bv lFuhr. Martin fe Wic sc (198 &j). As solar 
refere nce we adopted A(Fe)Q— 7.5 dex ( Grevesse fe Sauvaj 



1 19981 ). The error estimate has been performed similarly to 
what we did for Li. Taking into account the typical uncer- 
tainty on the atmospherical parameters we estimate a sys- 
tematic error of about ±0.10 and ±0.15 dex for SGB and 
RGB stars, respectively. Random error obtained from the 
Cayrel formula and the continuum level uncertainty ranges 
from ±0.14 (for SGB) to ±0.03 dex (for RGB stars). 



4 EVOLUTION OF A(LI) 

The right panel of Fig. [T] shows the behaviour of A(Li) as 
a function of luminosity for the observed targets. Luminosi- 
ties have been derived from the best-fit BaSTI isochrone 
discussed above and according with the position of the star 
along the evolutionary sequences. The two horizontal dashed 
lines denote the magnitude and luminosity levels where two 
abrupt drops of A(Li) occur. Fig.[2]shows the values of A(Li) 
as a function of T e //. We can describe the evolutionary be- 
haviour of the Li abundance by identifying four regimes in 
Fig. [T]and d 

• TO (T eff ^5900 K; lg(L/L Q ) <0.5): stars located at 
the TO level have an average <A(Li)>= 2.30±0.02 dex 



(cr= 0.10). To check whether the distribution of Li abun- 
dances is consistent - given the observational random errors 
- with a uniform value A(Li)=2.30, we performed the fol- 
lowing test (see, e.g., an application of the same test to the 
analysis of initial H e abundances in GCs a n d of GC ages by 
Cassisi et ai~l 120031 ; IChabover et all 1 19961 ; ISalaris fe Weissl 



20021 ). For each individual star we have calculated a set 



van der Waals damping constants by iBarklem et al.l 
(2000) when available. The source code is available at 
http://wwwuscr.oat.ts.astro.it/castclli/sourccs/width9.html 



of synthetic Li abundances by randomly generating - us- 
ing a Monte Carlo procedure - 10000 abundance values, ac- 
cording to a Gaussian distribution with mean value equal 
to ^4(Li)=2.30, and a equal to the empirical random abun- 
dance errors. As random errors we have considered the errors 
due to the EW measurements (as described in Sect. 3.3) and 
the random errors in the T e // determination (as described 
in Sect. 3.1). This is repeated for all 35 TO stars, and the fi- 
nal 350000 synthetic values are combined to produce an "ex- 
pected" distribution for the entire sample, in the assumption 
that the detected Li abundance dispersion is not intrinsic, 
but due to the individual errors, assumed Gaussian. The 
statistical F-test was then applied to determine if this "ex- 
pected" distribution (well sampled by the extremely large 
number of elements), displays a variance that is statistically 
consistent with the observed distribution of 35 objects. As 
customary, we state that an intrinsic A(Li) range of values 
does exist if the probability P that the two distributions 
have different variance is larger than 95%. 

We find a value of P below 90%, and therefore we can 
conclude that the observed dispersion is formally consistent 
with the random errors associated to the Li estimates of each 
star, and that our measurements provide a uniform Li abun- 
dance A(Li)=2.30±0.02 for all objects with T eff ^5900 K 

This result points toward a general homogeneity of the Li 
abundance in M4 stars, in contrast to the large star-to-star 
scatter observed in other clusters (as 47 Tuc and NGC 6752) 
that indicates an intrinsic dispersion in A (Li). 

• SGB (5200 K < T eff <5900 K; O.5<lg(L/L ) <0.6): 
these stars are located along the SGB, from the TO to the 
base of the RGB. We found a decrease of A(Li) for increasing 
T e ff, with A(Li) ranging from 2.0 (in the hot SGB stars) to 
0.9 dex (at the end of the SGB). 
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Figure 2. Behaviour of A(Li) as a function of T e f f for individual 
stars. The error bars include both random and systematic errors. 
The dott ed line denotes the WMAP+BBNS value of A(Li) com- 
puted bv lCvburt et al] j200Sf) and the grey area the corresponding 
uncertainty. 



• lower-RGB (4800 K < T eff <5200 K; 
0.9<lg(L/Z/o) <1.7): stars in this luminosity range, 
corresponding to the portion of the RGB fainter than the 
RGB Bump (that occurs for this age and metallicity at 
lg(L/Z/0) = 1.73), share the same abundance of Li, defining 
a plateau <A(Li)>= 0.92+0.01 dex (<r= 0.08). 

• upper-RGB (T eff <4800 K — lg(L/L Q ) >1.7): these 
stars are brighter than the RGB Bump and display a sudden 
drop of A(Li), reaching A(Li)~0.23 dex. 

We remark t hat M4 is the sec ond GC studied so far, 
after NGC 6397 (|Lind et all 12009 bl ), where A (Li) has been 
measured along the TO-SGB-RGB sequence: in both clus- 
ters two drops of A (Li) have been clearly detec ted, consistent 
with the pattern already seen in Halo stars (|Gratton et al.l 
l2000h . These drops can be associated with two distinct mix- 
ing episodes: 

(i) the canonical mixing episode called First Dredge- Up that 
occurs when the convective envelope reaches inner regions 
where Li has been destroyed through the nuclear reaction 
7 Li(p,a) 4 He; this produces a drop of A(Li) by a factor ~ 15 
along the SGB; 

(ii) the extra-mixing episode after the RGB Bump, that 
occurs when the outward moving H-burning shell reaches 
the discontinuity in the H-profile left over by the con- 
vective envelope at its max imum penetration (see e.g. 
ICharbonnel fc Lagarde Il2010l ). T he stars at V~13.5 (where 
the RGB Bump is located, see e.g. iFerraro et al.|[l999h share 
the same abundance of the fainter RGB stars, while the 3 
RGB stars above this luminosity display a drastic decrease of 
A(Li). When the mean molecular weight barrier is removed 
after the RGB Bump, the surviving surface Li is destroyed 
and the stars exhibit A(Li) approaching zero. 

The Li abundance measured in M4 TO stars, 
A(Li)=2.30 dex, is in agreement with what found in 
all other GCs investigated so far (Fig. [3}. The metal- 



Figure 3. A(Li) mean values for TO stars in the GCs stud- 
ied so far, as a function of the initial cluster iron abundance 
measured from the RGB stars. Error bars indicate the disper- 
sion around the mean, normalized to the root mean square of 
the number of used stars. The dark grey area displays the range 
of value s obtained for Halo star s by different recent investi- 
gations |Bonifacio fc Molarolll997l; ICharbonnel fc Primas |2005| ; 



lAsplund et al.ll2006l : iBonifacio et alj|2007l : I Aoki et al, 1 120091 The 

solid line displays the theoretical prediction for the value of A(Li) 
in TO stars for an age of 11-12 Gyr, in case atomic diffusion 
is fully efficient, starting from an initial WMAP+BBNS value 
A(Li)=2.72. The vertical dashed line shows the boundary for the 
radiative levitation to start affecting appreciably the photospheric 
chemical abundances of Li (see Sect. 8 for details). The primordial 
A(Li) provided by WMAP+BBNS calculations and the associated 
uncertainty is also displayed (light grey region). 

poor ([Fe/H] = — 2.3 dex) GC M92 shows a mean value 
of A(Li)=2.36±0.05 dex, a= 0.18, (|Bonifaciol [2001 ). For 
NGC 6397 ([Fe/H] = -2.0 dex) di fferent stud i es poin t to- 
ward A (Li) bet ween 2.25±0.01 (jLind et all feoogbl ) and 
2.30+0.01 dex |Gonzalez Hernandez et al] 120091 ) in stars 
along the MS or around the TO. The value of A(Li) in 9 
TO star s of NGC 6752 ([Fe/H] =-1.4 dex) has been mea- 
sured bv lPasauini et al] (|2005l) who found A(Li)=2.24±0.05 
dex (<T=0.15). Fina lly, the metal-rich cluster 47 Tuc 
(|D'Orazi et al] l2010h shows a similar mean abundance, 
but with a large star-to-star scatter, A(Li)=2.26±0.14 dex 
(cr= 1.38). 

Even if small offsets can be present between these Li 
determinations (because of the different temperature scales 
and adopted NLTE corrections), all these measurements 
are consistent each other and with the Spite Plateau de- 
tected in the Halo dwarf stars (grey region) . Hence, we con- 
firm for M4 the Li content discrepancy with respect to the 
WMAP+BBNS results, already detected in the field stars 
and in other GCs. 



5 EVOLUTION OF [FE/H] 

Fig. [4] shows the estimated [Fe/H] as a function of T e // 
for the individual stars in our sample. The iron abundance 
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does not show any systematic trends with the temperature. 
By taking into account the entire sample, the average Fe 
abundance for M4 turns out to be [Fe/H] = — 1.10 ± 0.01 dex 
(cr=0.07). This val ue basically agrees with previous spec- 
troscopic analyses. Ilvans et all (|l999l ) derived an a verage 
iron content [Fe/H]= -1.18±0.01 d ex. iMarino et al.l (|2008[) 
obtained [Fe/H]= -1.07±0.01 dex, ICarretta et all (|2009d ) 
found a value [Fe/H]= -1.18±0.02 dex. Taking into account 
the adopted solar reference values, the absolute differences 
between our Fe determinations and the previous ones are of 
+0.06, -0.01 and +0.04 dex, respectively. By dividing the 
sample in 3 groups, e.g. considering separately 35 TO stars, 
16 SGB stars and 36 RGB stars, the derived average Fe 
abundances are [Fe/H]=-1.08±0.01 (<r=0.07), -1.13+0.03 
(o-=0.07) and -1.11+0.01 (a=0.07), respectively. Hence, no 
evident trend is detected between dwarf and giant stars. A 
F-test analogous to the one described for the TO sample has 
been performed on the [Fe/H] values of our whole sample of 
M4 stars, and has provided - as expected - no statistically 
significant indication of an intrinsic spread. 

Note that this analysis does not include any correction 
for the departures from LTE, due to the lack of NLTE cor- 
rections for individual Fe I lines at the metallicity of M4. 
Moreover, the precise order of magnitude of NLTE correc- 
tions for Fe lines is not well constrained, because some as- 
sumptions must be introduced (for instance the factor Sh 
needed to correct the H I collis ional cross-sections derived 
by the formula bv lDrawinl (|l968T ) and the accuracy and com- 
pleteness of the adopted model Fe atom). Although the de- 
partures from LTE are more critical at lower metalficities, we 
have t ried to use the corrections computed bv lGratton et all 
(1 19991 ) for a representative high excitation Fe I transition 0, 
for an overall metallicity [M/H]=-1.0 dex. These NLTE- 
LTE corrections are positive and quite small (~0.05 dex for 
dwarfs and ~0.03 dex for giant stars) and, if applied, they 
do not change the observed trend and are not capable to in- 
troduce a systematic offset in [Fe/H] between TO/SGB and 
RGB stars. 



6 SANITY CHECKS 
6.1 Median spectra 

To reduce the error sources and strengthen our result about 
the SGB stars (that have the lowest SNR), we averaged their 
spectra according to their individual T e ff. We obtained five 
spectra with enhanced SNR (~150 around the Li line) and 
repeated the entire analysis for each of them: normalization 
of the spectrum, re-computation of Teff by fitting H a wings 
and measurements of EWs for Li and Fe. The resulting new 
set of T e ff well resembles the average of the single spec- 
tra values within ~40 K, confirming the reliability of the 
individual Teff. Also, A(Li) and [Fe/H] measured in these 
median spectra agree within few hundredths of dex with the 
average Li and Fe abundances computed from the individual 
stars. 



6 Note that the majority of the adopted Fe I lines in our analysis 
has excitation potential > 3 eV. 
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Figure 4. Behaviour of [Fe/H] as a function of T e ff for the in- 
dividual stars in our sample. Errorbars include random and sys- 
tematic errors. 



6.2 The impact of different T e ff scales 

The Li abundance shows a plateau for temperatures larger 
than T e /f= 5900 K (see FigEJ) . The reliability of this trend 
is crucial, because stellar evolutionary models including pure 
atomic diffusion predict an increase of A(Li) before the di- 
lution due to the Firs t Dredge-Up (see e.g. the discussion in 
ISalaris fe Weisdl200ll ). 

The most relevant source of error in the A(Li) deter- 
mination is the assumption of the Teff scale. We therefore 
repeated the analysis by adopting two different photometric 
T e j j scales, derive d by means of the most recent (V — 7)p- 
T e ff relation s by iGonzalez Hernandez fe Bonifacio! (|2009l . 
GHB09) and ICasagrande et all (|2010l , C10), both based on 
the InfraRed Flux Method (IRFM). We remark that these 
T e ff scales are independent of the T e ff derived by the H a 
line and affected by different source of errors. In particular, 
T e ff inferred by IRFM are largely insensitive to uncertain- 
ties in model atmospheres (i.e. the adopted line broaden- 
ing recipe, the assumption of 1-dimensional geometry) at 
variance with those derived from the H a line, but they are 
affected by uncertainties on the photometry, absolute and 
differential reddening, photometric calibration and so on. 
We have derived (V — I)o colours for our stars from the WFI 
photometry, corrected for differential reddening effects (see 
Sect. 3.2). By adopting the GHB09 calibration, we find an 



average difference T^ B09 - 



T 



-77 K (<r= 125 K). Sim- 
ilarly, the temperatures computed by the C10 calibration 
provide a mean difference T°ff-T" f a f = - 29 K (tr= 123 K). 
These mean differences are small in absolute terms and also 
much smaller than the dispersion around their mean values. 
We therefore conclude that no significant difference is found 
between the Teff scales obtained from the Balmer line and 
from the IRFM. This finding reduces the major disagree- 
ment between previous generations of theoretical (V — I)o~ 
T e f f transformati ons and IRFM calibra tions (see the exten- 
sive discussion bv lWeiss fe Salarislll999l ) . 

We have however re-analyzed our stars by employing 
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also these two T e // scales, re-computing both gravities and 
microturbulent velocities. The adoption of IRFM T e f / scales 
does not change the global behaviour of our abundances, 
both Li and Fe. In particular, the iron content remains con- 
sistent between TO and RGB stars, without systematic off- 
set between the different group of stars. Variations of about 
~-0.05 dex or less are obtained for the average A(Li) of TO 
stars. 

These checks further reinforce the conclusion that the 
observed behaviour of A(Li) and [Fe/H] as a function of T e // 
is not an artifact due to an incorrect T e // scale. 



6.3 Comparison with lD'Orazi &c Marinol (|2010h 

The only other sample of Li abundance in M4 stars available 
to date concerns 104 RGB stars discussed by DM10. There 
is a small overlap between the 2 samples, with only 5 RGB 
stars in common. The agreement with their atmospherical 
parameters is satisfactory, despite different methodologies 
have been applied (i.e., for T e ff we use the comparison with 
the BaSTI isochrone, while DM10 relies on the excitation 
equilibrium method). Our T e /f values are higher by +36 K 
(a— 22 K) on average, and also the gravities show very small 
differences (of the order of 0.2 dex). The mean difference be- 
tween the iron abundances of these 5 targets is of only —0.04 
dex (a— 0.06 dex), including also a small offset in the so- 
lar zero-point. Although atmospherical parameters and iron 
abundances are fully consistent between the two studies, a 
large (~0.4 dex) difference is revealed in the Li content of 
3 stars, while for the other 2 objects DM10 provide A(Li) 
close to zero, lower than our measures. The discrepancies are 
not attributable to the quoted measurement errors, nor to 
the difference in T e f / (which leads to a difference in A(Li) of 
~0.02 dex only). Fig. [5] shows the Li line in our spectra for 
the five stars in common, compared with synthetic spectra 
computed with our abundances (solid grey curve) and with 
those by DM10 (dotted grey curve). Note that for the 2 stars 
with A(Li)~0 in DM10, the differences are consistent with 
the different SNR of the spectra, because in their spectra 
(with typical SNR 50-100) such weak Li lines are hidden in 
the noise envelope. For the other 3 stars the discrepancy is 
more critical. 

DM10 have employed ATLAS9 model atmospheres 
computed with th e option of approximate overshoo ting. 
As pointed out by Molaro. Bonifacio fe Primasl (|l995l ) and 
iBonifacio et al.l (|2009T ), the temperature stratification of AT- 
LAS9 model atmospheres with and without approximate 
overshooting can be very different in the optical depth range 
— 1< lg(i") <1, where the bulk of the weak/moderated lines 
is formed (including the Li line adopted in this work). We 
computed some ATLAS9 model atmospheres with approxi- 
mate overshooting and compared the synthetic Li line pro- 
file in the two cases. At the metallicity of M4, the differ- 
ence between the thermal structures of the two models is 
less severe than at lower metallicities (as those investigated 
by IBonifacio et al.ll2009h . and the adoption of the ATLAS9 
models with approximate overshooting leads to an increase 
of the Li abundance smaller than 0.1 dex. 

As a sanity check, we apply our procedure to the GI- 
RAFFE dataset of NGC 6397 stars studied by iLind et all 
(|2009bh . We retrieved the reduced data from the GIRAFFE 
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Figure 5. Portions of the GIRAFFE spectra (solid black) around 
the Li line for 5 RGB stars in common with DM10. Solid grey 
lines show synthetic spectra computed with our Li abundances, 
and dotted grey lines display spectra computed with Li abun- 
dances from DM10. The bottom-right panel shows the Li line 
in the spe ctrum of star 1044 7 belonging to NGC 6397, from the 
dataset bv lLind et al . (2009a). compared to a synthetic spectrum 
computed with their atmospherical parameters and Li abundance 
A(Li)=1.04 (see text for details). 



Archive maintained at the Paris Observatory and derived 
A(Li) for a sub-sample of their targets by using our proce- 
dure and their atmospherical parameters and EWs. Differ- 
ences of a few hundr e dths of dex with respect to the derived 
A(Li) bv lLind et al] |2009bT ) have been found. The bottom- 
right panel of Fig|5] shows a RGB star from the sample by 
ILind et "ah! (|2009bf )Tin comparison with a synthetic spectra 
computed with their atmospherical parameters and A(Li). 
From this test, we exclude any relevant systematic offset 
in our procedure to infer the Li abundance. As a further 
ch eck, we redetermine d the EWs of the Li line in a number 
of ILind et all l|2009bl ) stars, finding negligible differences. 
This confirms the reliability of our method to measure the 
EWs. These two tests give further support to the correctness 
of our A(Li) values for the RGB stars, but the discrepancy 
with DM10 remains unexplained and we postpone this issue 
to future investigations. 

Finally, we note that if our sample was affected by a 
systematic underestimate of ~0.3-0.4 dex for both RGB and 
TO stars, we would obtain A(Li)~2.7 dex at the TO of M4, 
a value much higher than the A(Li) measurements along the 
Spite Plateau and in general higher than all the Li measure- 
ments in Pop II stars of the last three decades. 



http:/ /giraffe-archive. obspm.fr/ 
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7 LI ABUNDANCE IN FIRST AND SECOND 
GENERATION STARS 

As already mentioned in Section 1, the current interpreta- 
tion of the anticorrelations observed in the GCs, involves a 
first generation of stars with efficient CNO, NeNa, and MgAl 
cycles in their interiors, whose chemically processed matter 
is transported to the surface, and ejected in the existing 
interstellar medium through mass loss processes. A second 
generation of stars is then formed out of this polluted mat- 
ter, so that nowadays individual GCs are populated by a 
combination of first and second generation stars. The na- 
ture of this polluting first generation stars is still debated, 
with two competing candidate s: intermediate-mass asymp - 
totic giant branch (AGB) star s dVentura fc D'Antona1[2 009) 
or fast rotating massive stars (|Decressin et alj|2007l ). In both 
cases Li is expected to be depleted, since the NeNa cycle oc- 
curs at temperatures higher than that of Li burning. How- 
ever, if AGB stars are responsible for intracluster pollution, 
they may have als o produced Li, through t he Cameron- 
Fowler mechanism (I Cameron fc Fowler 1 1 19711 ). so that the 
surface A (Li) may be able to increase to values similar to 
the initial abundances. Thus, it is crucial to compare the Li 
abundance in first and second generation stars, bearing in 
mind that only the former should be born from a gas with 
the cosmological A(Li). To this purpose, we study the be- 
haviour of A(Li) as a function of [O /Fe]. The [O / Fe] ab un- 
dance ratios have been derived by lLovisi et alj (|2010T ) by 
measuring the permitted O 1 triplet at 7771-74 A. Unfortu- 
nately, our dataset does not include Na transiti ons, usually 
adopted to well separate the two populations fsee lLind et al.l 
l2009bl : ILvQrazi et all [20101 . hereafter DM10), because the 
choi ce of the setups h as been done for different purposes 
f see lLovisi et al.ll2010h . 

Fig. [S]shows the distribution of TO and RGB stars (ex- 
cluding the stars located along the two Li drops, see Fig. [TJ 
in the [0/Fe]-A(Li) plane. For both dwarf and giant stars, 
the Li abundance does not show an y correlation with [O /Fe], 
similar to what found in 47 Tuc fp'Orazi et al.ll2010l ) and 
in M4 by DM10. Measuring the Spearman rank coefficient, 
we find Cs~ 0.05 and 0.06 for dwarf and giant stars, respec- 
tively. Such a test confirms the lack of correlation between 
Li and O content. 

Hence, the Li abundance appears to be very similar 
among first and second generation stars. Adopting as bound- 
ary between the two populations the median value of the 
oxygen distribution ([O/Fe]~0.30 dex), we find very similar 
average A(Li) values, both for dwarf and giant stars. This 
finding points out that Li production occurs in this cluster. 

In fact, if the progenitors of second generation stars 
are massive objects, they have destroyed their original 
Li content in their envelopes, whereas if AGB stars 
are responsible for intracluster pollution, they may have 
produced Li, through t he Cameron-Fowler mechanism 
(| Cameron fc Fowler 11971b . so that the surface A(Li) may be 
able to increase to values similar to the initial abundances. 
On the other hand, a point to recall is that the Cameron- 
Fowler mechanism to produce Li needs to some fine-tuning. 
The entire envelope of the AGB stars where Li has been 
produced must be ejected before the newly formed Li is de- 
stroyed. The uniformity of the A(Li) values around the TO 
and along the lower RGB points somewhat to Li produc- 
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Figure 6. Behaviour of A(Li) as a function of [O/Fe] for the 
stars of M4, excluding those located along the two Li drops (see 
Fig.CE}. 



tion in the polluters, as already discussed by DM10. In fact, 
since our sample of stars at a given T e j j contains a random 
mixture of stars belonging to both first and second genera- 
tions, if the polluters destroy Li we would expect an intrinsic 
spread in A (Li). Conversely, if the polluters produce Li, a 
quite fine tuned scenario is required to achieve a uniform Li 
abundance for both first and second generations of stars in 
M4. In any case, several uncertainties (as the mass-loss rates 
and the Li yields) dramatically affect these kind of models. 



8 THEORETICAL MODELS 

To compare our measured Li and Fe abundances with 
the results from theoretical stellar evolution models - 
with and without including the effect of atomic diffu- 
sion - w e have employed t he st ellar evolution code de- 
scribed in iPietrinferni et al.l (|2006t ) . More in detail, we have 
considered an initial chemical composition with Y= 0.25, 
[F e/H] = —1.1, and t he sa me q-enhanced meta l mixt ure as 
in lPietrinferni et al.l (|2006t l and lSalaris fc Weiss! (|200ll ). with 
an average [a/Fe]= +0.4 dex and with initial A(Li)= 2.72, 
corresponding to the WMAP+BBNS Li prediction. We have 
then computed the evolution of models with masses of the 
order of 0.85-0.89 M GJ and TO age equal to 12 Gyr (for 
models without diffusion) and 11 Gyr (for models with dif- 
fusion). These objects are assumed to be the counterpart 
of stars currently evolving along the SGB and RGB of the 
cluster. The evolutionary sequences were computed starting 
from the pre-MS phase, until they reached T e // ~4700 K 
along the RGB, i.e. below the lower T e ff limit of our ob- 
served stars. 

Diffusion was included by solving the Burg- 
ers equations for a multicomponent fluid following 
iThoul. Bahcall. fc Loebl (|l994 ). and we considered explic- 
itly the evolution of H, He, Li, C, N, O, and Fe abundances. 
The other elements have been assumed to diffuse like Fe. 
Our calculations account for the effects of temperature, 
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gravity and molecular weight gradients, while radiative 
levitation is not included. 

As for the computation of the radiative opacity, given 
that for this metallicity all heavy elements are approxi- 
mately equally diffused, the abundance ratios are not signif- 
icantly affected, and the effects on the stellar opacity are ac- 
counted by interpolating among a-enhanced opacity tables 
of different Z. Superadiabatic convect ion is treated accord - 
ing to the mixing length formalism of I Cox fe Giuli I (|l968h . 
The value of the mixing le ngth has been chosen fro m a solar 
model calibration (see e.g. IPietrinferni et al.|[2006l ). Finally, 
the rates for the Li de struction reactions are taken from the 
NACRE compilation (|Angulo et al.ll 19991 ') . 

We remark that the only element transport mechanisms 
accounted for in our models are convection (with boundaries 
fixed by the Schwarzschild criterion and no overshooting) 
and atomic diffusion. No additional mixing process able to 
operate in RGB stars brighter than the RGB bump, or tur- 
bulent mixing that moderate the effect of atomic diffusion 
are included. 

Given that our analysis focuses on the theoretical inter- 
pretation of the observed Li and Fe abundances, it is impor- 
tant to assess whether neglecting the radiative levitation - 
that tends to counteract the effect of temperature, gravity 
and molecular weight gradients - can affect our conclusions. 
Regarding the pred icted surface Li ab undances, Fig. 3 o f 
IRichard et all (120021 ) and Figs 3 and 4 of lKorn et ail i|2006br ) 
show that for metallicities above [Fe/H]~ —2 the effect of 
neglecting radiative levitation is negligible. 

As for the surface Fe abundances, the inclusion of ra- 
diative levitation at metallicities above [Fc/H]~ —1.3 causes 
only small changes that do not affect any of the conclu- 
sions of our study. To give a quantitative example, we cal- 
culated a [Fe/H] = — 1.61, 12 Gyr isochrone, and compared 
the surface [Fe/H] values with t he co rresponding predic- 
tions in Fig. 9 of IRichard et al] (|2002T ). that include also 
the effect of radiative levitation. We find a depletion of 



[Fe/H] close to the TO (T e 



6300 K ) of .35 dex, 



compared to ~0.30 dex from 
Around the TO (T eff ~ 6200 K) of a [Fe/H] 



Richard et all (|2002l ') results. 

1.31, 12 Gyr 



isochrone we find a [Fe/H] depletion b y 0.27 dex, compared 
to ~0.25 dex from the corresponding IRichard et al.l (|2002l ) 
isochrone. 



9 DISCUSSION 

This section presents a detailed comparison between our Li 
abundance measurements and theoretical models, with the 
aim to assess the consistency between the initial A(Li) for 
M4 stars and the WMAP+BBNS value. 



9.1 Comparison with the theory: the role of the 
diffusion 

From the discussion in Sect. 4 it is clear that also M4 ex- 
hibits the discrepancy between the A(Li) value measured at 
the TO and the WMAP+BBNS estimate. We first investi- 
gate the role played by the atomic diffusion, that has been 
suggested by other authors as a po ssible solution to th e prob- 
lem, at least in case of NGC 6397 (|Korn et alj|2006ah . Fig. [7] 
and [8] display the comparison between the measured surface 
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Figure 7. Behaviour of A(Li) as a function of T e j j for our sam- 
ple of stars grouped in eleven T e f f bins. The curves indicate the 
theoretical predictions with (dashed line) and without (solid line) 
diffusion. Dotted lin e denote the primordial A(Li) computed by 
ICvburt et al.l (j2008Ti and the grey area the corresponding uncer- 
tainty. 



Li and Fe abundances (grouped in 100-150 K wide temper- 
ature bins only for sake of clarity [f|) and the predictions 
of evolutionary models with and without atomic diffusion, 
from the TO to the lowest temperature limit of our sample. 
For sake of comparison, we also plot the WMAP+BBNS 
A(Li) value used in the initial composition of our models, 
and the related uncertainty. The value A(Li) at the TO in 
the calculations without diffusion is very similar to the ini- 
tial value (lower by only 0.03 dex) because the pre-MS Li 
depletion at these metallicities and for the TO mass of M4 
is negligible. This value of A(Li) stays constant (within the 
uncertainties) from the TO to a point along the SGB with 
T e jf ~5900-6000 K; below this temperature, the surface Li 
abundance drops steadily, for the inward extending base of 
the convective envelope is reaching increasingly hotter lay- 
ers where Li has been burned. This surface depletion ends 
when the convective envelope attains its maximum depth at 
the beginning of the RGB. 

In the model with diffusion, the surface Li around the 
TO is d epleted by about 0.2- 0.25 dex, in agreement with re- 
sults bv lRichard et al.l (|2002l . see their Fig. 9). With decreas- 
ing temperature along the SGB, the surface A(Li) starts to 
rise again due to deepening of the convective region, that re- 
turns to the surface some of the Li diffused outside the enve- 
lope along the MS. A(Li) reaches a peak when T e // ~5800- 
5900 K, with an abundance within ~0.1 dex from the initial 
one. Moving to lower effective temperatures, the surface Li 
abundance starts to drop, because the convective envelope 
reaches regions where Li had been burned due to the higher 



8 Note that the adoption of different binning steps does not 
change dramatically the results. 
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Figure 8. Behaviour of [Fc/H] as a function of T e ff for our sam- 
ple of stars grouped in eleven T e ff bins (same binning of Fig.[7]l. 
The curves indicate the theoretical predictions with (dashed line) 
and without (solid line) diffusion. 

temperatures of the stellar matter. As for the case without 
diffusion, surface depletion ends when the convective enve- 
lope attains its maximum inward extension. 

Fig. [9] helps to better clarify this behaviour. The Li 
abundance profile inside a model representing the cluster TO 
stars is displayed as a function of the temperature, for cal- 
culations with and without diffusion. The model calculated 
without diffusion displays a constant A(Li) even beyond the 
convective envelope boundary, located at log(T)~6.02, down 
to layers where the temperature is high enough to burn Li. 
The very small (~ 0.03 dex) difference with respect to the 
initial value is due to Li-burning during the pre-MS stage. 
When the model evolves along the SGB, convection deep- 
ens. As long as the bottom of the convective envelope stays 
within the temperature range where the Li profile is con- 
stant, the surface abundance does not change. When the 
base of the convective envelope reaches log(T)=6.3-6.35, it 
engulfs Li-depleted layers and the surface A(Li) decreases 
steadily, until the envelope reaches its maximum extension 
at the end of the SGB phase. At this point the surface de- 
pletion ceases and the surface A (Li) predicted by the models 
stays constant during the following RGB evolution. In case 
of models with diffusion, below the base of the convective 
envelope A(Li) first increases, due to the accumulation of Li 
diffused below the convection zone during the MS. When T 
reaches the Li-burning temperature, A(Li) drops fast, as in 
the model without diffusion. Along the SGB evolution, the 
base of the convective region engulfs first the layers with 
increased A(Li), and this causes the observed surface abun- 
dance increase after the TO. When convection goes deeper 
than log(T)=6.3-6.35, the behaviour of the surface A(Li) be- 
comes similar to the case without diffusion. It is important 
to notice that the observed constant value of A(Li) along 
the RGB - before the RGB bump mixing episode that is 



not accounted for by the models - provides an additional 
constraint on the Li content at the TO, that is rarely con- 
sidered in the majority of the chemical analyses. In fact, 
given that the maximum depth of the convective envelope 
attained at the base of the RGB is independent of the exact 
value of the initial Li abundance, the predicted A(Li) along 
the lower RGB is determined by the internal Li abundance 
profile at the TO. 



As for the trend of [Fe/H] with T e //, the model 
without diffusion displays essentially a constant [Fe/H] 
along the whole evolution, equal to the initial value (the 
effect of the First Dredge-Up, that decreases slightly the 
surface H-abundance, is negligible). Instead the model 
with diffusion has an abundance lower by 0.21 dex around 
the TO, due to the diffusion of Fe below the convective 
envelope during the MS phase. Along the SGB [Fe/H] 
increases steadily, for the increasing depth of the surface 
convective region brings back to the surface most of the 
Fe previously diffused from the convective envelope. When 
the model reaches the RGB, the surface Fe abundance is 
only very slightly lower than the initial value, and stays 
constant thereafter. We estimate that the T e ff of the TO 
stars should be increased by at least 300 K in order to 
lower the estimated iron abundance to values consistent 
with the predictions of models including atomic diffusion, 
while keeping the same T e // for the RGB stars. However, 
the resulting T e // for the TO stars would be incompatible 
with the T e ff of theoretical isochrones, with or without 
diffusion. 



Some considerations can be drawn from the comparison 
between the observations and the theoretical models: 

(1) we measure a uniform value of [Fe/H] along all the evolu- 
tionary sequences, that is reproduced by the model without 
diffusion. Taken at face value, this suggests that a process 
able to inhibit atomic diffusion from the convective envelope 
must be at work in the M4 stars; 

(2) both the models with and without diffusion, that start 
from a WMAP+BBNS initial Li abundance, predict A(Li) 
for TO stars higher than our measurements in M4 (see 
Fig.0; 

(3) there is no observational evidence within the current ob- 
servational errors for the A(Li) peak before the SGB drop, 
predicted by the models with atomic diffusion (see Fig. [TJl ; 

(4) both the theoretical curves in Fig. [7] match reasonably 
well the value of T e ff (~5900 K) where Li dilution due to 
the penetration of the convective envelope starts, and the 
temperature (T eff ~5200-5300 K) where the First Dredge- 
Up ends. However, a clear offset in A(Li) is recognized. In 
fact, the models predict A(Li)~0.30 dex higher than that 
observed among the SGB and RGB stars. This offset can- 
not be erased by adopting different T e f / scales (s e Sect. 6.2). 
It is in teresting to note that also in NGC 6397 (|Korn et al.l 
l2006al . see their Fig. 1), the Li abundance among the RGB 
stars is lower than that predicted by the models tuned to 
fit the TO and SGB data (R05) by ~0.2 dex (despite this 
vertical offset has been not discussed by the authors). 
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Figure 9. Li abundance profile (in terms of the difference with 
respect to the initial value adopted in the calculations) as a func- 
tion of temperature, for a model corresponding to M4 TO stars, 
including (dashed line) and neglecting (solid lino) diffusion. The 
dot-dashed line denotes the abundance profile needed to repro- 
duce the measurements of A(Li) in the cluster (see Fig. HOI . 



9.2 Comparison with the theory: the role of the 
turbulence 

A quite debated route to solve the discrepancy between the 
WMAP+BBNS Li abundance and the measurements along 
the Spite Plateau is to assume that, besides atomic diffusion, 
a turbulent mixing process occurring below the convective 
region is at work, partially or totally counteracting the ef- 
fect of diffusion from the envelope. The models developed by 
R05 account for atomic diffusion (including radiative levita- 
tion) plus an ad iioc parametrization of turbulent mixing 
that limits the settling of Li (and other metals), and at the 
same time does not transport matter into the convective 
envelope where Li has been burned by nuclear reactions. 
In particular, the mode ls computed with r eference tempera- 
ture log(T)=6.O0 allow iKorn et all l|2006af ) to reproduce the 
observed TO and SGB A(Li) measurements in NGC 6397, 
starting from a primordial A(Li)=2.54 (lower than more re- 
cent estimates from WMAP+BBNS calculations) and the 
small [Fe/H] difference (~0.17 dex) they measure between 
TO and RGB stars. 

Let's consider the internal Li profile for the model cal- 
culated with diffusion, as that displayed in Fig. [5] Any tur- 
bulence that brings back into the convective envelope some 
of the Li accumulated around log(T)=6.3 will increase the 
surface A(Li) and, at the same time, lower the Li abundance 

9 This parametrization makes use of a turbulent diffusion coeffi- 
cient that is 400 times larger than the diffusion coefficient of He 
at a reference temperature T, and varies with density (p) as p~ 3 . 



peak below the envelope. As a consequence, also the early 
increase of surface A(Li) during the SGB phase tends to be 
erased. This mechanism works well for NGC 6397, where 
A(Li) for TO stars is larger than the value predicted by 
atomic diffusion models, and the early increas e along the 
SGB i s also lower than expected from diffusion (|Korn et alj 
l2006al ) . The same efficiency of turbulence (i.e. the same ref- 
erence temperature log(T)=6.0 in the parametrization by 
R05) however, cannot work for M4 and probably other clus- 
ters. Fig. [3] highlights the problem. The solid line repre- 
sents the value of A(Li) in TO stars, as a function of the 
metallicity, predicted by our models including atomic dif- 
fusion (without radiative levitation) for a reference age of 
11-12 Gyr, and starting from the WMAP+BBNS value of 
A(Li). As discussed in Sect. 8, the lowest metallicity for 
which these Li predictions are adequate is [Fe/H]~ —2.0 dex, 
since below this value radiative levitati on affects appreciabl y 
the predicted surface Li abundances (|Richard et alj|2002h . 
When TO stars are located above the solid line (as in the 
case of NGC 6397), atomic diffusion during the MS is too 
efficient in depleting the Li surface abundance and one can 
reproduce the observations only by including some form of 
additional turbulence that brings back to the surface some 
of the Li diffu s ed out of the convective envelopes, as done 
bv lKorn et al.1 (l2006an . For objects located below the solid 
line - like M4 - however, the inclusion of turbulent mixing 
with the same efficiency as in NGC 6397 worsens the dis- 
crepancy, because it increases the surface A(Li) predicted 
by fully efficient diffusion. 

This suggests that - within the assumption of an initial 
A(Li) equal to the cosmological value - for the three highest 
metallicity clusters in Fig. [3] including M4, the turbulence 
has to reach deeper regions where Li burning occurs (i.e. one 
needs a higher reference temperature in R05 parametrization 
of turbulence), and Li-depleted material is dragged to the 
surface, thus lowering A(Li) below the value expected from 
atomic diffusion. The same turbulence could also account for 
the constant Fe abundance observed in our M4 star sample, 
feeding back to the envelope iron that diffuses below the 
convection zone during the MS phase. 

The observed trend of A (Li) with T e ff constrains the 
shape of the Li stratification in M4 TO stars. The value of 
A(Li) measured at the TO and its constancy (within the 
measurement errors) along the early SGB (without the Li 
peak typical of the Li profile in case of diffusion), plus the 
constraint given by the efficiency of the Li-burning reactions 
require a profile like the one displayed as a dot-dashed line 
in Fig. [5] We have calculated the SGB and RGB evolution 
starting from a TO model with this 'observationally con- 
strained' Li profile and the abundances of the other elements 
as in the calculations without diffusion, in the assumption 
that turbulence has mixed back into the fully mixed enve- 
lope the amount diffused below the convection boundary. 
The resulting surface Li evolution (Fig. 1 1 Op reproduces very 
well the abundance pattern up to the RGB bump, including 
the A(Li) values along the early RGB. 

9.3 Can we really detect the signature of 
diffusion? 

As discussed in Sect. 4, the observed distribution of A(Li) 
in TO stars is formally compatible with a constant Li value, 
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within the measurement errors. In the previous analysis we 
have therefore considered a constant A(Li) for all TO stars. 
The question we address here is the following: given the ran- 
dom errors on our estimates of A(Li), is it possible to detect 
the trend of A(Li) with T e // predicted by models with fully 
efficient diffusion (dashed line in Fig. [7} ? To this purpose, 
we have created a synthetic sample of A(Li) measurements, 
drawing randomly 35 values of T e ff from the Li-isochrone 
displayed in Fig. using a flat probability distribution for 
T e ff 1? 5900 K. Gaussian errors equal to our average ran- 
dom observational errors for the TO stars have been then 
applied to this synthetic sample, and the a F-test as de- 
scribed in Sect. 4 has been performed. This procedure has 
been repeated 100 times, and for each of these 100 samples 
the F-test returned a probability below 90% that a spread in 
A(Li) is present. This means that our measurements of A(Li) 
cannot rule out the existence of a trend consistent with pre- 
dictions from models with fully efficient diffusion. If this were 
to be the case, the observations are matched for an initial 
A (Li) ~0.20-0.25 dex lower than that the WMAP+BBNS 
value. 

However, a similar test performed on the [Fe/H] values pre- 
dicted by the models with diffusion for whole T e // range 
of our data, shows that the F-test should be able to detect 
the expected spread of values, which is contrary to the result 
for our observed sample. We have verified that the predicted 
[Fe/H] depletion has to be reduced by at least ~70% in or- 
der to become undetectable by the F-test, given our random 
measurement errors. This means that the efficiency of dif- 
fusion has to be strongly suppressed by come counteracting 
turbulence and that the curvature in the Li abundance as 
a function of T e f / for TO stars has to be much lower than 
what predicted by fully efficient diffusion, if not completely 
vanishing. 

9.4 The initial Li content of M4 stars 

The discussion in the previous subsections and the fit 
displayed in Fig. [TO] do not yet allow to identify univocally 
the pristine A(Li) in M4, because this value is strongly 
dependent on the details of the adopted turbulent mixing, 
that seems to play an essential part in the interpretation of 
the spectroscopic abundances. Two extreme situations can 
be envisaged to obtain the model plotted in Fig. 1101 

(i) a pristine A(Li)= 2.72 dex and a turbulent mixing able 
to reach inner regions where Li is burned, thus dragging 
to the surface Li-poor material. This solution requires a 
higher reference temperature in R05 parame trization of 
turbule nce with respect to that proposed by iKorn et al.l 
|2006al lbl), that does not reach stellar regions where Li is 
depleted. This solution allows to solve the discrepancy with 
WMAP+BBNS. 

(ii) a pristine A(Li)~2.35 dex and a turbulent mixing 
qualitatively similar to that of R05 with log(T)= 6.09, 
showed in their Fig. 7. This efficiency would erase the Li 
peak at log(T)= 6.3 (see Fig. [9]), by increasing the envelope 
A(Li) back towards the initial value. In this case another 
possibility to explain the discrepancy with the cosmological 
Li abundance - apart from inadequacies of the BBNS 
standard model that could lead to a n incorrect estimate of 
the c osmological A(Li) value (see e.g. IJedamzik fc Pospelovl 
l2009h - is to invoke that an amount of primordial Li 
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Figure 10. Behaviour of A(Li) as a function of T e f f for the stars 
of M4 (grouped in T e f f bins as in Fig.0, in comparison with the 
theoretical model computed by adopting the Li profile showed in 
Fig. [9] as dot-dashed line. 

could have been processed by Pop ulation HI stars, before 
Population II stars were formed (|Piau et all |2Q06j). The 
implication in this case is that Population II stars are 
not ideal to infer the cosmological Li abundance. While 
this could partially (up to 0.3 dex) solve the discrepancy 
between the primordial A (Li) and the initial abundance 
in M4, it also requires a complex and fine tuned scenario. 
Even if Population III stars have depleted the Li abundance 
before the cluster formation, some form of turbulent mixing 
has still to be active below the convective envelope of M4 
stars, to erase the A(Li) peak around log(T)=6.3, that 
otherwise would produce an increase in the surface A(Li) 
at the beginning of the SGB evolution. 

In light o f these considerations, the solution proposed 
bv lKorn et al.l (|2006al ) , based on R05 parametrization of tur- 
bulence with reference temperature log(T)=6.0 and a initial 
Li abundance of 2.54 dex, may not be a universal solution 
for the discrepancy with WMAP+BBNS calcu lations. This 
value of log(T) adopted bv lKorn" et al.l (|2006ah can account 
for the observations of NGC 6397, but would make the dis- 
crepancy worse for M4. Instead, in case of M4 (and possibly 
also of NGC 6752 and 47 Tuc) one needs a turbulent mix- 
ing much deeper than in NGC 6397, that burns a sizable 
fraction of the initial Li. 



10 CONCLUSIONS 

We have measured the Li and Fe content for a large sample 
of stars along different evolutionary phases (TO, SGB and 
RGB) in the GC M4. The main conclusion of this work can 
be summary zed as follows: 
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• The Li content in the TO stars of M4 well resembles 
the Li abundance already observed in other GCs and in 
the Halo stars, thus confirming the discrepancy with the 
WMAP+BBNS predictions. 

• Qualitatively, the general behaviour of both Li and 
Fe with T e ff is reproduced by models without diffusion, 
starting from the WMAP+BBNS A(Li) value but we 
recognize a systematic offset of ~+0.3 dex between the 
theoretical and the observed Li content along the SGB and 
RGB stars. 

• The behaviour of A(Li) as a function of T e // in 
TO/SGB stars is incompatible with the theoretical pre- 
dictions of the models including atomic diffusion, starting 
from the WMAP+BBNS predicted value; This suggests 
that an additional turbulent mixing below the convective 
zone needs to be included in the models. 

• The general homogeneity and the value of A(Li) in 
the TO, early SGB and early RGB stars constrains to 
the shape of the Li stratification in M4 TO stars. A good 
match between theoretical models and the global behaviour 
of the Li abundance as a function of T e / / can be obtained 
by assuming the WMAP+BBNS predicted value plus 
diffusion and a turbulent mixing reaching regions where 
Li is burned. The discrepancy with the WMAP+BBNS Li 
abundance cannot b e solved by invoking t he same reference 
temperature for the iRichard et al. | 120021) pa r ametri zation 
of turbulence, as proposed by iKorn et al ] l|2006al lbh to 
explain the Li abundance measured in NGC 6397, because 
in that case the turbulent mixing does not reach layers 
where Li is burned. 

• Different (lower than the WMAP+BBNS predicted 
value) pristine A(Li) values for M4 are possible, depending 
on the adopted reference temperature for the turbulent 
mixing. This uncertainty is mainly due to our poor under- 
standing of the true nature of the turbulent mixing that 
counteracts the atomic diffusion. 



The study of the Li content in M4 points out that 
the Li problem is far from being solved and new efforts 
are needed, both to enlarge the number of studied GCs 
and to refine theoretical stellar and cosmological models. 
Spectroscopic measurements of A(Li) below the TO would 
be crucial to put further constraints on the primordial Li 
in GC stars, given that moving down along the MS - even 
in case of fully efficient diffusion - one expects to find Li 
abundances progressively closer to the primordial value. 
Moreover, a deeper understanding of the nature of the 
turbulent mixing occurring outside the convective region is 
needed to disentangle this problem. 
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